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This series of lectures consists of two parts. In the first part the foundations of 
Q\ ' perturbativc and non perturbative formulation are discussed. The ambiguity in 

' the definition of vacuum condensates is then analyzed. In the second part the 

, symmetry patterns of the deconfining phase transition are discussed. 

7— I 1 

' 1 Introduction 

<D • 

' I shall preliminary review the basic formulation of quantum field theory, to 

^\ , clarify what is meant by perturbative and non perturbative. This will also 

help to understand the foundations and the limitations of lattice QCD, which 
is the main existing non perturbative approach. 

■ 1-1 Basic Feynman path. 

■ A quantum system is described by a set of canonical variables q(t), p(t) and 
y— i , by the Hamiltonian H. q(t),p(t) is a short symbol for q(i, t), p(i, t) the index 

' i running over the degrees of freedom. When the index i is not specified what 
we say will apply to any system, independent of the number of degrees of 

+^ ■ freedom. In a quantum field theory q(i, t) are the fields 

'T 1 ; q(i,t) = (p a (x,t) (1) 

^ '" the index being i = (a,x). 

, The system has a continuous infinity of degrees of freedom. 

Solving a field theory means constructing a Hilbert space on which q, p act 
as operators, obeying the canonical commutation relations and the equations 
of motion. A ground state must exist to ensure stability. 

The usual approach consists in splitting the Hamiltonian in a free Hamil- 
tonian Hq and an interaction H[ 

H = H + H! (2) 

Hq is quantized and solved exactly, and the resulting Hilbert space is assumed 
to describe also the interacting system. Hj is an operator in that space and is 
treated as a perturbation. This approach is called perturbative quantization. 
A more fundamental formulation can be given in terms of Feynman path 



- 1 — i 

X 



integrafl. 
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The basic quantity is the action 

S[q] = 



Cdt 



(3) 



with C the Lagrangean. S' is a functional of the fields. The equations of 
motion are the stationary points of the functional 

ss 



5q l {t) ' 

The knowledge of the partition function 







Z = 



l[dq(i,t) 



exp[iS(q)] 



or better of the generating functional 



Z[J]= J \{dq{i,t) exp{i[S(q) + Y, J dtJi(t) qi (t)]} 



(4) 



(5) 



(G) 



provides the solution of the theory. 

Indeed the connected time ordered correlators of the fields in the ground 
state, are expressed in terms of Z[J] as 

6 n Z[J] 



•\T(q(i 1 ,t 1 )...q(i n ,t n ))\0) 



— (-i) n 

Z SJ(i 1 ,t 1 ) . . .SJ(i n ,t n ) 



(7) 



In terms of the correlators the Hilbert space is determined, and the solution 
of the theory is known (reconstruction theorem). 

Z is a functional integral, with a continuous infinity of integration vari- 
ables. By definition it is dcfined as the limit of a sequences of ordinary inte- 
grals, computed on a discretized space time, on sets of points which become 
dense in space time in the limit. _ 

The above statements are clearly illustrated by the following argumentEI. 
Consider the transition amplitudes 



(q',t'\ q ,t) 



(8) 



between eigenstates of the g's at different times. Any physical amplitude will 
be known if all the amplitudes (||) are known, the q's being a complete set of 
commuting operators. 

Eq.(||) can be rewritten 

{q'\e- iH V~%) (9) 
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Putting S = 

e -iH(t'-t) = [ e -^<5]" +1 

as n — > oo, J — > 0. 

By inserting complete sets of states, we get 

(q',t'\q,t) = ]T (Q'\e- iH5 \qn)(qn\e- iHS ...e- iHS \ qi )( qi \e- iHS \ q ) 

9n— 9l 

2 

We shall consider for the sake of simplicity H's of the form H = v + V(q), 
but what follows can be easily adapted to momentum dependent interactions. 
Then, by use of the Baker-Haussdorf formula 

e -iHS _ e -i£-s e -iV{ q ) ^ + ^2^ ^ 



(q k+1 \e- iHS \q k ) 



^ 



(11) 



and 



(q't'\qt) = lim / TT c% cxp 

J fc=l 



= lim 



f| <*7fc e 



iS n (q) 



(12) 



fc=i 



Eq.'s ( pl|JT^ ) allow to continue analitically to imaginary times (Euclidean 
space time) by the Wick rotation xq = it. The limit of the Euclidean ampli- 
tude when x' — > T, x — > —T, and T — > oo can be evaluated as follows 

(g , T|9-T) B = J] (g'|e- iH |^n 1 )(£;„ 1 |e- [S(9)+ - / - 9] | J En 2 )<^ 2 |e-? ff | g ) 

E ni ,E n3 

= e-* £ " 1 fe„ 1 (g')e-* £ " 2 V' £ „ 2 (<z)(^ 1 |e-[ s ^ +J ^]|i?„ 2 } (13) 

E ni ,E„ 2 

The external current has been taken to be zero outside the time interval 
( — 1'' l')' anc ^ the are ^h e wave functions of the states |£"„) which are 
eigenstates of energy. As T — ► oo only the ground state survives and, apart 
from irrelevant multiplicative constants which cancel in eq. (^) , eq. (|l3|) defines 
the euclidean version of Z [J], and selects the true ground state of the theory. 
The Lagrangean defines a field theory if the limit 8 — ► (n — ► oo) exists 
(ultraviolet limit), and if the limit T —> oo exists (infrared limit). 

Note that after Wick rotation Z is in fact the partition function of a 
statistical system. The infrared limit exists if the theory ha a mass gap, or, 
in the language of statistical mechanics if the thermodynamical limit exists. 
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1.2 Lattice QCD 

Lattice formulation of qcd! is a clever construction of approximants in the 
sequence of discrete integrals which define the partition function. 

The theory is defined on a square lattice in space-time, and the building 
block is the link 

U^n) = e iagoA ^ n) (/i = 1...4) (14) 

which is the parallel transport from the site n of the lattice to the neighbor- 
ing site in the direction fj,. The action is expressed in terms of the parallel 
transport along an elementary square, in the plane /iz/, the so called plaquette 

S=J2^eTr\l-^uA (/3 = ^r) (15) 
In the limit of zero lattice spacing it tends to the continuum action 

S a ~,--J2 G % G U ( 16 ) 

A critical point exists at (3 — > oo (g 2 — > 0) where the (3 function of the 
theory tends to zero by negative values (asymptotic freedom)0, and like in 
any statistical model QCD is defined as a field theory when approaching that 
point. 

At sufhciently large /3, by renormalization group arguments the lattice 
spacing a, measured in physical units, becomes smaller and smaller 

a((3) ~ eM-boP) (17) 

— bo being the first coefficient in the perturbative expansion of the (3 function. 

At sufhciently large (3 the correlation length becomes large with respect 
to lattice spacing, and the coarse structure of the lattice irrelevant: the limit 
of eq.([l2|) exists. 

On the other hand if the lattice is kept large compared to the correla- 
tion length also the thermodynamical limit is reached, since the theory has a 
fundamental scale built in. 

This argument indicates that QC'D exists as a field theory, even if a formal 
proof of that statement does not yet exist in the sense of constructive field 
theory. 

At sufficiently large (3, and large enough lattice a good approximant is 
obtained for Z, in the sense of the limit ([l2|), and in particular the true 
ground state is selected, as shown in eq.(fl3l). 
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The formulation is gauge invariant and needs no gauge fixing and no 
Faddeev-Popov ghosts. 

It would be nice if analytic solutions could be given of the lattice regu- 
larized theory. Unfortunately nobody has been able up to now to produce 
thcm. Numerical simulations are the usual technique to compute physics 
from lattice. Euclidcan QCD is a statistical mechanics, and like in statisti- 
cal mechanics Montecarlo techniques can be used to compute the regularized 
Feynman integrals and correlation functions. 



1.3 Perturbative expansion. 

The action S can be split in two parts: a term So which is bilinear in the 
fields and a term Si which includes the rest 

S = S + St (18) 

Sq can be written as 

a0 

where a — i, t spans the degrees of freedom and time. Then exp[— S] is 
expanded in powers of Si as 

cxp[-S] = cxp[-5 ] Y, L -T L ^ 2 °) 



n'. 

n=0 



and 



Z = jY[dq a exp[-5] = f [] dq a e- S ° f] { —^- (21) 

a a n—O 



S j is in general a polynomial of the fields and by inserting eq.(|2C 
(T(q ai . . . q a J = — Y[dq p exp{--q a V^ qi3 )q ai ...q an Y^ — ~j— ( 22 ) 

J p n=0 

Ali terms of the expansion are gaussian integrals of the form 
1 f 1 

7? / H ^ eXP (~ ^IcPallP)^! ■'■q<Tn ( 23 ) 

which can be computed and give 

Y, V atiat V at3 ^V a ^_ iatn (24) 

pairs 
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where the sum runs over all possible pairings of indices. Eq.(|2J) is nothing 
but Wick's theorem with V the propagator. 

In QCD such scheme of quantization produces an S matrix describing 
scattering of quarks and gluons, and the ground state is Fock vacuum. The 
thcory can be renormalized producing for any observable quantity a power 
series of the renormalized coupling constant a s , with finite coefficients 

O = £0 P a£ (25) 

However the series is not convergent, not even as asymptotic series. This 
reflects an instability, probably due to the inadequacy of Fock vacuum as 
ground state. It is phsically obvious that the theory describing scattering 
of particles which are H$¥er observed cannot be a good starting point for a 
perturbative expansionQ'cl. 

It is an empirical and not fully understood fact, however, that a few terms 
of the expansion (|2^) give a good description of phenomena involving short 
distances. 



1.4 Instanton vacuum. 

An attempt to go beyond perturbation theory was developedi after the dis- 
covery of classical solutions of QCD with finite action (instantons), having 
non trivial topologyu. 

In the computation of an ordinary integral the saddle point approximation 
implics a sum over all the stationary points of the phasc: 

f dx e*™ ~ W d(Sx n )e^ +I ^ 5 < ^T" ( 26 ) 

J x„ J k' 

A/„ = / - f(x n ) - ^-5x1 

Also in the computation of the Feynman integral by the saddle point approx- 
imation all the stationary points of the action should be considered. 

The usual perturbative expansion consists in retaining only the trivial 
stationary point q = 0, and the saddle point expansion is nothing but eq.(|2^). 
If other minima exist, with finite action, then, the gaussian approximation 
reads 



So«9-V ^ (A/ ra ) fc 
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S [q n ] is the action of the classical configuration q n , Sq n = q — q n , AS 

S-S°[qJ i-^ 



2 5q a 5qp 



The existence of minima, i. e. of stable classical solutions. is due to topol- 



ogy. Finite action requires that, in the euclidean space tim 

FuvFuv — O ( —r 



lons. 
imeB, 



i. e. that the gauge field on the 3d sphere at infinity S3 is a pure gauge. 
This defines a mapping of S3 into the gauge group. To each continuous con- 
figuration an integer is associated which indicates how many times the group 
manifold is swept in this mapping. If dU is the volume element of the group 
normalized to 1, (/ dU = 1), and da 11 the volume element of S3 



dU = K^da* 1 / K^da* 1 = n (28) 
By Gauss theorem, if d^K^ — Q{x) 

' Q(x)d i x = n (29) 

Q(x) is known as topological charge density. 

Explicit form of the instanton was produced in ref.Q. 

At the end of the 70's a big effort was done to include non perturbative 
effects via eq.(p7|), i. e. by taking into account extra terms beyond the trivial 
stationary point q = 0. The hope was that small instantons would dominate 
and hence a few terms of the perturbative expansion around them would be 
a good description of hadronic physicstl. 

The attempt failed due to difficultics in controlling the infrared problems, 
and to the difficulties in explaining confincmcnt. It is known, however, that 
instantons play an important rolc in understanding chiral properties of the 
theorya. 

1.5 Non perturbative approaches. 

We have no real knowledge of non perturbative methods in field theory. Only 
numerical simulations on lattice allow to extract information on the mecha- 
nisms of the theory, e. g. on confinement. 

A taily non perturbative idea which has proved to be fruitful is the N c — > 
00 limitEI. 

The idea is that the main properties of gauge theories are independent 
of N c , and the same as for N c — > 00 (1/N C —> 0), with g 2 N c = X fixed. 
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Corrections for finite iV c , are perturbations which do not alter the structure 
of the theory. 

Fermion loops are of order Nf /N^ 2 1 {N f is the number of light quarks, 
n the number of vertices). Except for n — 2, the graphs which enters in the 
defmition of the (3 function, all the rest should be a small correction. Lattice 
simulation do indeed indicate that, except for a rescaling of the masses due 
to the /3 function, numerical determination in quenched approximation differ 
in general by few percent (~ 10%) from the determinations in full QCD. 

The 1/N C limit gives an explanation of the so called U(l) problem, relating 
it to a quantity, the topological susceptibility of the vacuum x>L2l. Lattice 
determination of x confirm the validity of the approachtll. The same idea 
proves a useful guiding principle in understanding confinement (see sect.3). 



2 Non perturbative effects in QCD. Vacuum condensates. 

2.1 S. V. Z. Sum rules. 

Consider the Wilson operator product expansion of the correlator of two cur- 
rents 

W(x) = T (r(x)f(0)) = C?{x) I + Cr^G a a pG a a0 

+C^{x)m^ + ... (30) 

In principle this expansion is a theorem ip_perturbation theoryEl, but is pre- 
sumably valid in QCD at short distancesEj. In Fourier transform and renor- 
malizing 

U MV (g) = J i 1 !^ (0^(^(^(0)10) = 

= (q»qv ~ q 2 9^) [n(g 2 ) - n(0)] (31) 
n(q 2 )-n(0) = C 1 + %G 2 + %G v , + ... (32) 

Q 2 = (o | ^-^-G^G ^ 1 0) has dimension [m 4 ], is the vacuum expectation value 
of the dilatation anomaly and is known as gluon condensatc. 

Gij) = (O|mi/i'0|O) & l so nas dimension [m 4 ], is related to the breaking of 
the chiral symmetry and is known as the quark condensate. 

Ci, C^Cip are dimensionless constants, apart from log's, and depend on 
a s . G2, G^ are not defined in perturbation theory. A dispersion representation 
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for n( g 2 ) = r% 2 ) - n(o) is 



% 2 ) = V f ^^L_ (33) 

7 4m 2 s s~ q z + i£ 

with i?(s) the well known ratio 

R(s) = °^-^ s ] (34) 

In perturbation theory, at each finite order only C\ gets contributions and is 
a constant modulo logs. a ~ l/s, modulo fogs. 

Physics, however, is not scale invariant: resonances break scale invariance 
by powers, and condensates should describe that breaking. This is the basic 
idea of SVZ sum rules where the equality "on average" is exploited 

-,.r*« , C ,4 C ,4 C , (35) 

J 4m 2 s s-q 2 +is g 4 q A 

Equality "on the average" means that the two sides are not equal at each 
value of q 2 , but suitable averages are equal, which emphasizc Ihe region where 
the two description are expected to be approximately valio£3. The C coef- 
ficients are computed to some order in perturbation theory. A successfull 
phenomenology results ouLof which Gi and (mipip) can be extracted. 
From a recent reviewt-3 we quote 

G 2 = (0.24 ± 0.011) GeV 4 (36) 
W)iGeV^-0.13GeV 3 (37) 

Two questions arise 

a) Is the Wilson O.P.E. eq.© well defined? 

b) Can the constants G2, C,p be determined by a non perturbative formu- 
lation of the theory, like lattice? 

The answer tp the first question is that it is ambiguous, due to the presence of 
renormalonsH. The answer to the second is yes, modulo the same ambiguity. 

In spite of this ambiguity, however, keeping the leading terms of the per- 
turbative expansion, sum rules do work, and so does the lattice determination, 
which agrees with phenomenology. 

We will sketeh the argument for the first question, we will show that it 
applies in the same way to lattice determinations, and will present anyhow 
the successful results from lattice, which agree with the values ( |36| , j37|) . 
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The perturbative expansion of any observable, e. g. the correlator of two 
currents like II^„, eq.(32), is finite order by order in perturbation theory, i. e. 



d(a a 



(38) 



where a s is the measured couplmg constant, and arc finite numerical 
coefficients. It can be shown thaua 



C[ n) oc n\n a b r ' 



(39) 



where a does not depend on the operator, b does depend on it. Of course the 
dependence ( |39| ) on n makes the series non convergent. However the series 
could be considered as the asymptotic expansion of a convergent series if the 
Borel transform 

C 1 (b) = J2c[ n+1)b - (40) 

n=0 



is convergent. In that case, since 



J b n e- b/a °db = a n s +1 n\ 



the asymptotic expansion 

dia.) -Ci(0) = / dbe- b ' a >d{b)db (41) 
Jo 

is just eq.(|3^). 

Most probably this is not the case for C\{a s ). 

The resummation of a subseries of graphs show^that G\ (a s ) is ambiguous 
in principle by terms which mimic the condensatecJ. 

Consider the QCD correction Hi(q 2 ) to the correlator, corresponding to 
the following graphs 
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Replacing a s by its running value at one loop a s (k), l/a s (k) = — + I 2 - ln 



niiq2) = 7?J Wf asG{k > q2) (42) 

nning value at one loop o 
amounts to resum a subset of graphs and gives 

At small fe's, G(fc, q 2 ) is finite and equal to AN c /q 2 and hence (^) has a 
divergence at -C /i. The divergent part is 

1 , 60 , k 

, . . . h — m — 

rf J \a s Z7r /i 



ni (g 2 )^9^/fc 3 dfc(- + ^ln-) (44) 



The prescription of the singularity is not defined. Accordingly there is an 
ambigui ty 

ni(g ) - — 2 — 4 e * 45 ) 

which mimics the terms proportional to the condensate. 

The argument is only an indication, since the subset of graphs considered 
is not even gauge invariant. If the indication is correct it means that the Borcl 
transform is singular, and condensates cannot be defined 

As already anticipated, in spite of that sum rules work and give a consis- 
tent phenomenology. 

-A^ similar problem arises when one tries to extract G2 from Iatticec30t3. 
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Consider the gauge invariant field strength correlatoiS 

V^ p< , = (0\Tr {G^{x)S{x, Q)G pa {0)S\x, 0)} |0) (46) 



where 



r 1 \ rpa r~ia 

a 

S c (x,0) = Pexp jf A fJ ,(x)dx> x J 



T a are the generators of the gauge group and Sc(x, 0) is the parallel transport 
from to x in the fundamental representation. S depends on the path C: our 
choice for C will be a straight line. 

'Dfjtv^v can be considered as the split point regulator of G2, a la Schwinger. 

The OPE gives 

2W ^ n ^(I)+C i G 2 + ... (47) 

x^0 X 

< G2 > is defined if the term oc -K can be unambigously subtracted. Again 
like for the correlators, higher order perturbative terms containing log's can 
sum up and simulate a constant term, like G2, with an intrinsic ambiguity in 
the prescription. 

The most general form of the correlator compatihle with Poincare invari- 
ance is expressed in terms of two independent scalarstj D(x 2 ), D\(x 2 ) 

T^fip^a = {g^gpc - g^crg vp ) [d{x 2 ) + Di(x 2 )] + (48) 

dD 1 



~t~ y&fA&ugpa •^p'^agvp % p% cr 9 pv '^U'^pOfiaJ 



dx 2 



Callingj—for simplicity the direction of Xp two other invariants can be 
definedliZl 

z— 1 



T>± (x 2 ) = - J2 = V + P i ( 5 °) 

i<j 

Dii and £>x can be computed on the lattice0. 

Their behaviour is shown in fig.l for quenched theory. 
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V{x 2 ) = G 2 e 



(51) 



allows to extract from the data G 2 and A. 
For quenched S U (3) 

G 2 = (0.15 ±0.03) GeV 4 A = (0.22 ± 0.02) fm (52) 

For full QCD and 4 quarks at different values of the quark masses (in units 
of inverse lattice spacing) the result is smaller than the quenched one by an 
order of magnitude and the correlation length is larger 

am = 0.01 G 2 = (0.015 ± 0.005) GeV 4 A = (0.34 ± 0.04) fm (53) 
am = 0.02 G 2 (0.031 ± 0.005) GeV 4 A = (0.29 ± 0.03) fm (54) 

By use of the relation0 

dG 2 24 



dm, 



11 



(55) 



and of the known values of (qq) an extrapolation can be made to the physical 
quark mass, givingtS 

G 2 = (0.022 ± 0.005) GeV 4 



which compares satisfactoraly^ith the sum rule determination, eq.(^ 
Old results for SU(2) giveEl 

G 2 = (0.10 ±0.04) GeV 4 



(56) 



(57) 
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A few words on the technique used in refffl. Correlators measured on 
a statistical ensemble of configurations are very noisy due to local quantum 
fluctuations-of the fields. A local cooling procedure is used to reduce quantum 
fluctuationsE 2 !. After a number n t of cooling steps the size A of the fluctuations 
which are cooled away follows a diffusion equation 

A 2 = Dn t (58) 

If the distance to explore is large enough (e. g. 3p41attice spacing) cooling can 
eliminate local noise and reduce the fluctutionsEZEa. 

Field correlators as determined from lattice are useful inputs in the sto- 
castic approacbo to QCD vacuum, for which I refer to the lectures of Y. 
Simonov of this school. 

An alternative procedure to determine G2 is to use ±lie plaquette as a 
regulator, and measure the average of the plaquette itsehO A perturbative 
contribution appears, which is presumably non Borel summable, which must 
be somehow subtracted to isolate the non perturbative part. Again there is 
ambiguity. However if the perturbative subtraction is made at some order, 
the resulting G2 fits what is found by other methods. .— . 

By a similar procedure the quark condensate can be extractedll3 from the 
correlator 

m(tp(x)S(x, 0)^(0)) ~ -4 + C 2 (m^) (59) 
The result is, at (3 = 5.35 



|(^)| 7 - (205 ± .21)MeV am = 0.02 ( — = 0.65(3)) (60) 

Tflr, 



l' 

(160±12)MeV am = 0.01 —=0.57(2) (61) 

\m p ) 

This approach is like sum rules: it works but it is not understood why 

3 Confinement of colour. 

3.1 Confinement. 

The particles corresponding to the fundamental fields of the thcory, quarks 
and gluons, have never been observed in nature (confinement of colour). 

In nature the relative abundance of quarks to nueleons has an upper 
bound 

^ < 10^ 27 (62) 
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coming from Millikan like analysis of ~ 1 gr of mattet- 
The expectation in the absence of confinement isE3 

^ > 1(T 12 (63) 
n p 

a factor 10 15 larger. 

A factor 10 -15 cannot be explained by fine tuning of a small parameter: 
the only natural explanation is in terms of symmetry. We shall adopt this 
prejudice as a starting point. 

Lattice provides evidence that confinement is present in QCD. The well 
known area law of Wilson loopsEj, W(R, T) ~ exp(— RTa) defines the static 
potential between q q as 

V (R) = Ra (64) 

a is known as string tension: a ^ indicates confinement, meaning that an 
infinite amount of energy is necessary to puli apart q and q. 

Finite temperature QCD can be obtained from the Euclidean partition 
function, with a time extension l/T, and periodic boundary conditions for 
bosons, antiperiodic for fermions. 

On the lattice this is done by using a lattice N^Nt of extension N$ 3> Nt 
in each space direction and 

N T a = I (65) 

with a the lattice spacing. 

By renormalization group arguments at large enough /?, a depends on /3 

as 

a ( / 3)^_Lexp(-6 /3) (66) 
Al 

— bo is the first coefficient in the expansion of the /3 function, I/Al the physical 
length scale of the theory. 

Because of asymptotic freedom a decreases exponentially to as (3 — > oo. 

The temperature is given by eq.(|6^j6^) as 

T = exp(6 /3) (67) 

and is an increasing function of (3. At low temperature a ^ 0. By incrcasing 
T (or (3) a decreases and at £QWi£ temperature T c it goes to zero. T c is called 
the deconfining temperaturea' 
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An alternative signal for the transition is given by the Polyakov line 

V = (T I exp(iA (n, x )dx )) (68) 

which is the parallel transport along the time axis from uq to infinity and 
back via periodic boundary conditions. 

It can be shown that P is rclatcd to the chemical potential of an isolatcd 
quark, /i g as 

V = exp(-^) (69) 

In the confined phase [i q diverges and V = 0, while V ^ in the deconfined 
phase. Such behaviour is indeed observed on. the lattice, again providing 
evidence that there is a deconfining transitioncj. 

V is called an order parameter, being ^ in the weak coupling regime. 
a is called a disorder parameter, being ^ in the strong coupling regime. 

A drawback of a and P as order parameters is that, in the presence 
of quarks they lose their meaning. Indeed the string should break at some 
distance creating qq pairs, and a is not defined. As for V, it is an element 
of the centre of the group Zj, and signals the breaking of Z% symmetry in 
the confined phase. In the presence of quarks Z% is not a symmetry. The 
order parameter usually assumed in full QCD is the chiral symmetry breaking 
order parameter {ipip}. A priori, however, chiral symmetry is not related to 
confinement, although arguments show that it is if the philosophy of N c — > oo 
is correct. In the same philosophy, however one should be able to definc 
an order (or disorder) parameter which is independent of the presence of 
quarks, at least for not a too large number of light quarks, which would spoil 
asymptotic freedom. 

Because of asymptotic freedom the confined phase corresponds to disor- 
der, the deconfined phase to order. In the language of statistical mechanies 
one should then look for the symmetry of the system dual to QCDcZl. 

3.2 Order disorder duality. 

Order disorder duality is a property of d + 1 dimensional systems, with d- 
dimensional extended configurations fi having non trivial topology, and car- 
rying a conserved topological charge. 

In the weak coupling regime the system is usually deseribed in terms of 
local fields and the symmetry of the ground state by order parameters ($), 
which are vacuum expectation values of some field A phase transition can 
oecur at higher couplings, where ($) — > and the system becomes disordered. 
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However a dual description can be given of the system, in which the 
configurations fj, become local fields, the original $'s become extended objects, 
the coupling g — > l/g, and hence the disordered phase looks ordered, with 
order parameter (/i), (/i), the order parameter of the dual phase is callcd a 
disorder parameter (which means order parameter of the dual). _ 

The prototype example of duality is the 1 + 1 dimensional Ising modelEa. 
The model is described by a field a which takes values ±1 on the sites of a 
2d square lattice, and is exactly solvable in the thermodynamical limit. The 
interaction is 

"EE^M^) £ = exp|/3^5>W cr ( i + ^} (70) 

k=l i l fc=l i J 

The model has 1 dimensional configurations with nontrivial topology, the 
kinks, the topology being identified by the values of a at x — ±oo at fixed t. 
A kink is a highly non local configuration with 

a = —1 x < xq a = +1 x > xq (71) 

A dual variable a* can be defined which has also values ±1 on the sites of the 
dual lattice, whose sites corresponds to the links of the original one. 

At fixed t a* is — 1 if the link joins neighbourghs sites on which a has the 
same sign, +1 otherwise. 

In terms of a* a link looks as a* = +1 at xq, a* = —1 everywhere else. 
The kink becomes a local excitation in the language of a* . 

In ref.cS it was shown that the partition function in terms of a* is again 
an Ising model with the correspondcncc 

Z[a 7 (3} = Z[a*,P*} (72) 

and 

Shlh2/r = simW (73) 

The mapping to the dual sends (3 — > 1/(3 i. e. the weak coupling region into 
strong coupling region and viceversa. The high temperature phase which looks 
disordered, is ordered and (p*) is its order parameter. The order parameter 
of the dual phase is called a disorder parameter. 

Other known systems which admit a dual description and a disorder pa- 
rameter are 

a) The (2 + l)d XY model, (liquid He^). Here the field is 8 (i) an angle, or 
better —d^9{i). The topological excitations are 2d vortices. Topology 
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is given by their winding number and 

<j> A - dx 

which is an integei. if the field 6 is single valued. Vortices condense in the 



disordcrcd phasc 



l'.i 



b) The (2 + l)d Heisenberg modelEl Here the topological excitations which 
condense in the vacuum above the Curie temperature, when the mag- 
netization goes to zero, are non abelian vortices, the instantons of the 
2 dimensional 0(3) a model. Topology is the winding number of the 
mapping of the 2 dimensional sphere S2 on S O (i). 

c) In (3 + 1 \_dimensional gauge theories the topological excitations are 
monopolesEH, whose topology comes from the mapping S2 on SU(2)xlt 
Si on U(l). Examples are the (3+1) dimensional U(l) gauge theoryEJ, 
the (3 + 1) dimensionaLpon abelian gauge theory, the (3+1) dimensional 
supersymmetric QCDE2l. 

To understand disorder one can adopt two strategies. 

1) Perform the transfonnation to dual. This has been done in ref.ll for, the 
Ising model, in ref.E3 for the U(l) compaet Villain aetion, in refr3 for 
the SUSY QCD. 

2) Guess the symmetry of the dual phase and write the disorder parameter 
(/i)4n terms of the original local fields $. This line has bee: 
ref.E-3 on the continuum and by the Pisa group on the lattio 



'l •H' !• 



The advantage of strategy n. 2, which we shall adopt, is that the symmetry 
can be explored without going through the hardly manageable procedure of 
construeting the dual. 

The general principlc is to create by an operator /i, constructcd in terms 
of the local fields <I>, a dual excitation carrying non trivial topological charge, 
and investigate dual order by directly computing (fi), the disorder parameter. 

This is done by the field theoretical analog of translation operator, 

e lpa \x) = \x + a) (74) 

The operator0 



t) = CX P 



J d 3 xll k (x,t) $ k {x,y) 



1 

k 



(75) 
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with Hi(x, t) the conjugate momentum to <fi(x, t) and y) the field config- 
uration produced by a topological excitation located at y. In the Schrodinger 
representation, where the fields (p i (x, t) are diagonal, 

H{y,t)\<pi(x,t)) = \(pi{x,t) + ®i{x,y)) (76) 

A topological configuration $ has been added to the original one. 

Care is needed to adapt eq.'s( [75| , [76|) to compact theories, where the fields 
cannot be arbitrarily shifted. Another important point is that (/z) is only 



defined in the thermodynamical limit. Indeed by eq.(75), if the number of 
dcgrecs of frecdom is finite (/i) is an analytic function of /3, and cannot vanish 
identically in the ordered phase without vanishing identically. This is only 
possible when V — > oo and Lee- Yang singularities develop in the partition 
function. Ali that can be mastered. We have computed (/i) for a number of 
models (in the 3d XY model, the Heisenberg magnet, the (3 + l)d compact 
U (i)) where the symmetry of the dual is known. We have thus a tool to 
explore the deconfining transition in QCD. 

As we shall discuss in QCD the precise symmetry of the dual system is 
not known. Configurations candidate to condense in the confined phase are 
monopoles. It is indeed a long standing ideaHcj, that confinement could be 
produced by dual Meissner effect. Dual here is intended in the sense electric- 
magnetic, i. e. that the rolcs of electric and magnetic fields arc interchanged 
with respect to an ordinary superconduetor. The chromoelectric field aeting 
between a qq pair is expected, by this mechanism, to be confined by dual 
Meissner effect into Abrikosov fiux tubes, with energy proportional to the 
length. 

The string tension a is then the energy of the tube per unit length. 

The disorder parameter would then be the v. e. v. (/i) of any operator which 
creates a monopole. (/i) ^ signals that monopoles condense in the vacuum, 
or, under very general assumptions, that the system is a dual superconduetor. 

The preliminary question is then to identify the monopoles. 

3.3 Monopoles in non abelian gauge theories. 

Let us first consider SU (2) gauge group, for the sake of simplicity. Let $>(x)cr 
be any local operator in the adjoint representation, in the usual notations. 

We can define the unit veetor parallel to $ as $ = $/|$|, everywhere in 
a configuration, except at the zeros of $(a;). 

Let us defineciJ the field operator F^ v as 

F MV = $<3^ ~^-{D^AD„$)-$ (77) 
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F^v is a colour singlet, defined everywhere except at zeros of $(x), and is 
invariant under non singular gauge transformations. 
Here 

G»» = d„A u - d v A» + gA„ A A v (78) 

is the gauge field strength and Z? M = <9 M — igA^A is the covariant derivative. 

The coefficient in front of the second term is chosen in such a way that 
the bilinear terms in A^ cancel between the two terms. Indeed one has 

Ffu, = $ ■ {d^A v - dvA») - i (d M $ A dA) ■ $ (79) 

Definc 

^ = JW^ (80) 

and 

if = (si) 

Then j*^ is zero wherever is defined, i. e. everywhere except at the zeros of 
$, and is in any case conserved since F* is antisymmetric. 

W? = (82) 

Eq.([S2"|) identifies a conservation law, a magnetic U(l) symmetry. 

One can operate a gauge transformation bringing $ along a fbced axis. 
By eq.(|79|) then, in all nonsingular points 

F„ v = dpA* - d^Al (83) 

is an abelian field. That is why such transformation is called an abelian 
projection. 

It can be shown that at the singularities of $ a non zero magnetic cur- 
rent is present. The regular field is a magnetic field of a monopole sitting 
at the singularity. The singularity of the gauge transformation which per- 
forms the ahelian projection generates the Dirac string which ensures flux 
conservatiorci 

The abelian projection of the t'Hooft Polyakov classical monopole 
configurationcJ does indeed demonstrate that F^ v is the field of a static Dirac 
monopole. 

In general, for an arbitrary choice of a conserved magnetic charge can 
be defined. So a functional infinity of conserved monopole charges exists. 

Which of them condense in the vacuum to produce dual superconduetiv- 
ity, if any? 

There are in literature two different attitudes on the problem. 
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a) (t'Hooft)cZl Ali monopoles, whatever the choice of are physically equiv- 
alent. Ali of them should condense in the confined phase, to produce dual 
superconductivity of the vacuum, and go to normal state above T c . 

b) Some monopole are more equal than others. The maximal abelian pro- 
jection identifies the relevant monopolesEl 

We have checked the two alternatives on the lattice as explained in the fol- 
lowing sectionsSj. 

A similar construction works for SU (3), where two different monopoles 
species can be defined for each operator $ = r ^-$ a in the adjoint representa- 
tion, carrying two different magnetic chargescU. 

4 Some details on the disorder parameter (/i). 
The recipe of eq. (|7^) amounts to say that 

Indeed 

M = \ J tll W e M~PS - PAS} (85) 
Let us consider £7(1) for simplicity 

exp[-/3A5] =exp{-/3[(n 0j (%+& 4 )-n 0j (%]} p=\ (86) 

Here l/ebi(x — y) is the vector potential describing the field of a monopole. 
n i = [cos(#oi) — 1] is the plaquette 0, i. 

In the weak coupling 0Qi = Ei, the electric field and 

AS~-\[eEi}h 

gives exactly the definition ( f75| ) for /i. 

A similar definition can be given for QCIj23. 

Instead of computing (/i), which is affected by wild fluctuations, it proves 
convenient to compute 

P=^ ln <M> ( 87 ) 



From eq.(84) 



p=(S)s~(S + AS) s+AS (88) 
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In terms of p 



(fi) = exp 



(89) 



since by dcfinition (/j) = 1 . 

In the thcrmodynamical limit one expects in finite temperature QCD that, 
if monopole condcnsation is related to confinement, 



lim (jm) = 

V -^oo 

lim ( M ) ? 

V — »00 



(3>(3 C 
0<0c 



(p) ~ (/3 C - /3) 5 /3 ~ /3 C 



(90) 
(91) 
(92) 



At finite volume the qualitative behaviour of (/z) is shown in fig.2, the quali- 
tative behaviour of p(/3) in fig.3. 



< JI > 




Fig.2 
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▲ 




Fig.3 

For /3 > P c (weak coupling) p can be computed in perturbation getting, 
in terms of the volume 

(p) = cxp(- Cl y + c 2 ) ci > (93) 

(p) tends to zero as V — > oo. 

For f3 < (3 C p stays finite at — > oo, or (/i) is ^ 0. 

Around /3 C the phase transition takes place, which is known to be second 
order for SU '(2) gauge group, weak first order for SU (3). This means that 
the correlation length goes large for both, since both of them look like second 
order. On dimensional grounds 

M^M(f,~) t*(/3 e -/3)- v (94) 

£ being the correlation length and a the lattice spacing. 
Near the critical point a/£ ~ and 

M = /x(|,0) = /(i 1 / 1 '(/3 c -/3)) (95) 

Eq.(p4|) has been used. 
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This gives the scaling law for p = ln /3 



L 1 ' 



P=F((L^(p c -f3)) 



(96) 



From that scaling v and /? c can be determined, together with the critical index 
6. 

If condensation is related to confinement the determination should be 
consistent with other methods. 

The behaviour of p for monopoles defined by different choices of <I> is shown 
in fig.4: it looks the same for all of thcm, in agreement with the statement of 
t'Hooft that all monopoles are physically equivalent. 
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Also the two different monopoles corresponding to the same abelian pro- 
jection look identical. 
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Fig.5 
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The quality of scaling law for SU(2) is shown in fig.6, for the optimal 
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determination of [3 C and v. 
The result for v is 



v = 0.62 ± 0.05 (97) 

and for S 

S = 0.18 ±0.06 (98) 

For SU (3) the transition fig.7 
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is first order with 

;/ = 0.33 ±0.07 S = 0.54 ±0.04 (99) 

Monopoles do condense in the vacuum below T c , supporting thus the mcch- 
anism of confinement by dual superconductivity. Monopole identified by dif- 
ferent abelian projections behave in the same way. As a check the finite size 
scaling analysis of provides determination of /3 C and of the critical indices 
consistent with other methods. 

The disorder parameter (/i) can be defined and determined also in full 
QCD, contrary to a and P. Preliminary results givc evidence for monopole 
condensation also there. 
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The mcthod has been checked on 2d Ising model, on 3d XY model, on 
Heisenberg magnet, on U(l) with similar results, consistent with completcly 
diffcrent determinations. 

5 The symmetry of the confined phase in QCD. 

What we have learned from lattice is that QCD vacuum is a dual supercon- 
duetor in the eonfined phase, and goes normal at the deconfining transition. 
And this indcpcndcntly of the choice of the magnetic charge. 

This solves the problem that in a given abelian projection some colored 
states could be neutral with respect to the abelian electric charge, and hence 
not confined. 

However the situation is not elear. We do not understand yet the symme- 
try pattern of the ground state, except for the faet that any magnetic charge 
condensate is a good disorder parameter. There must exist for sure a more 
synthetic way of deseribing this situation. 

An attractive feature of our order parameter is that it is independent 
of the presence of quarks, and this fits the idea that the N c — > oo theory 
is a good basic deseription of physics, 0(1/N C ) correetions and higher being 
perturbations which do not alter its basic features. 
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